Devices with electrically tunable topological insulating phases 
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Solid-state topological insulating phases, characterized by spin-momentum locked edge modes, 
provide a powerful route for spin and charge manipulation in electronic devices. We propose to 
control charge and spin transport in the helical edge modes by electrically switching the topological 
insulating phase in a HgTe/CdTe double quantum well device. We introduce the concept of a 
topological field-effect-transistor and analyze possible applications to a spin battery, which also 
realizes a set up for an all-electrical investigation of the spin-polarization dynamics in metallic 
islands. 



The original prediction of the quantum spin Hall phase [l|, [2| has generated a renewed interest in topological 
phases in solid state systems. In the following years, the realization of the topological insulator (TI) phase has been 
theoretically described and experimentally observed in two-dimensional (2D) HgTe/CdTe quantum wells (QWs) [3|,|I|. 

ID helical edge modes are present at each boundary between a 2D TI and a normal insulator (NI), including the 
vacuum. Due to the topological protection against non-magnetic disorder [5|, [6| and the spin-filtered nature of such 
ID channels, TI edge modes are extremely promising for spintronics. A crucial point, however, is to devise a reliable 
method allowing to control transport in these helical channels. Several authors have proposed quantum interference 
and the Aharanov-Bohm effect as a possible route to control charge and spin transport in 2D TIs [TT-KUj. However, 
the most straightforward route would be to change the TI into a NI and completely turn off the helical channel. Two 
recent proposals show that an electrically tunable TI phase can indeed be obtained in an InAs/GaSb type-II QW [HI 
and in a HgTe/CdTe double quantum well (DQW) [l6|. 




FIG. 1: (Color online) (a) Isometric sketch of a HgTe/CdTe DQW device with a back gate and two distinct top gates (left and 
right). In the ON state, top gates induce a gate-bias domain leading to a TI/NI interface (channel) where helical edge modes 
are found. Source (S) and drain (D) leads, placed along the interface between L and R top gates, collect charges from the edge 
modes. The lateral surface of the DQW is specifically treated to ensure negligible edge transport, (b) Schematic description 
of a TI/NI interface for direct gate polarization with indication of the helical spin transport of edge states, (c) Reverse gate 
polarization leading to opposite spin transport of the channel. 



HgTe DQWs are driven from NI to TI by the application of a inter- well potential bias \V\ > Vc, where Vc is a 
critical value on the order of the gap of the individual QWs [l6j|. As shown in Fig. [Ha), a left (L) and right (R) top 
gate are employed to generate an inter- well potential bias domain (PBD), with V = Vl and V = Vr in the L and 
R region, respectively. The device is turned ON when the system realizes a TI/NI interface between the L and the 
R regions and helical edge modes run along the PBD line. Hence, we distinguish between a direct gate polarization 
accodingly to Vl > Vc > Vr and Fig. QJb), where spin up electrons run from source to drain, and a reverse gate 
polarization Vl < Vc < Vr with opposite spin transport properties [Fig. QJc)]. The device is turned OFF when the 
L and R regions belong to the same topological class, i.e. Vl,Vr < Vc (both NI) or Vl,Vr > Vc (both TI). Source 



* Electronic address: michetti@physik.uni-wuerzburg.de 



2 



and drain electrodes, which could for example be obtained by diffusing metallic atoms in the DQW system, lay close 
to the orifice between L and R top gates in order to collect electrons from the helical channel. 

The design in Fig. QJa) is studied to deal with a single helical channel (along the PBD), but helical edge states 
would also appear, if the DQW is in a TI phase, at the interface with the vacuum. We will assume the lateral surfaces 
of the DQW to be specifically treated in order to impede charge transport in such edge modes. Indeed numerical 
simulations on quantum spin Hall systems [17j suggest that transport of an helical edge mode at a physical edge can 
be suppressed by locally doping with magnetic impurities of random or in-plane magnetization. We will see that 
however such lateral edge modes, even when gapped out, play a relevant role in governing the spin properties of the 
system. 

The low-energy physics of a single HgTe/CdTe QW is known to be captured by the Bernevig-Hughes-Zhang (BHZ) 
model ||, which is, in first approximation, block-diagonal in the Kramers partner (spin) degree of freedom. For the 
spin up block the BHZ Hamiltonian reads Q 

written in the basis of the low energy subbands { +), |if 1, +) }. Conduction and valence bands of the TI for the 
spin up block are readily obtained by Eq. (pQ). Spin down eigenstates can then be obtained by applying the time reversal 
operator T = —iK, with K the complex conjugation operator. We adopt A = 375 meV nm, B — —1.120 eV nm 2 
and D = —730 meV nm 2 , which have been estimated by a comparison with the 8x8 Kane Hamiltonian fl8j |. and 
assume C = without loss of generality. The Dirac rest mass M depends on the QW thickness d. For d > 6.3 nm 
(M < 0) the system realizes a TI and the real space solution of Eq. (pQ) allows for helical edge states at any TI/NI 
boundary [3[. 

The BHZ model has also been generalized to the case of a HgTe DQW [Hj], where the two QWs are separated by a 
thin spacing layer of thickness t allowing for tunneling. The Hamiltonian for the spin up block has then the following 
structure 

(Kk) + \ h T (k) \ 

where the tunneling Hamiltonian hr(k) is given by 

hr(k)=( ,. Aj * + (3) 

v ' \a(k x - iky) A H i J w 

omitting terms of order 0(k 2 ). Numerical estimation of the tunneling parameters [l6| indicates that Aei ^> Ahi so 
that Ahi can be neglected. Moreover Aei and a can be taken constant as a function of k. Indicative estimates are 
reported in Table [J for t ranging from 4 to 7 nm. 
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Aei (meV) 


22 


12 


6 


3.4 


a (meV nm) 


18 


10 
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TABLE I: Tunneling parameters a and Aei for a HgTe/CdTe DQW, calculated at k x = k y = for several values of the spacing 
layer thickness t [161 ]. 



For defmiteness, we consider a structure of two equal QWs with M = 10 meV (QW thickness d « 5.5 nm), with a 
thin inter-well spacing layer of thickness t = 5 nm, ensuring a significant amount of tunneling, and consider a PBD 
where Vl = 30 meV and Vr = in the direct gate polarization [Fig.QJb)]. Fig.[2fa) shows the 2D bulk band structure 
of the system, where full lines refer to the L side and dashed lines to the R side. Note that the TI bandgap (L side) 
lays entirely in the NI bandgap (R side). The bandgap in the TI phase develops due to the inter-layer tunneling 
interaction and indeed, keeping fixed V = 30 meV, the TI gap is about Eg ~ 4, 2 and 1 meV for a spacing layer 
t = 5, 6, 7 nm, respectively. We solve now Eq. (j2j) in real space, imposing the continuity of the wave function and 
its normal derivative through the PBD line. When the device is ON, the system admits helical edge modes located 
at the TI/NI interface. The ID dispersion curve of such helical edge states is shown in Fig. Efb), together with the 
merging bulk bands. 

Let us now analyze ballistic transport in the helical edge states. Due to their spin- filtering nature, it is useful 
to define separate chemical potentials jii ajCr for the two spin species (a =t, i), with a = S, D for the source and 
drain contacts, respectively. Chemical potentials are kept in the DQW gap, so that at sufficiently low temperature 
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FIG. 2: (Color online) (a) 2D bulk dispersion curves for the left (L) and right (R) regions of the DQW device, which are in 
the TI and in the NI state, respectively. Bulk dispersion curves are degenerate in the spin degree of freedom, (b) Detail of 
the low energy spectrum with a plot of the ID helical edge modes at the TI/NI boundary (the PBD line). Edge states in full 
thin lines and in dashed thin lines belong to opposite Kramers blocks (spin). We employ the following parameters: t = 5 nm, 
M = 10 meV, V L 30 meV and Vr = 0. 



no charge is present in the bulk of the L and R regions. For simplicity, we consider charge accumulation in the edge 
states (which are considerably extended) to have negligible effects. The unidirectional charge current for the spin 
species a in the ID helical channel is given by 

e f Ea e 
4 C) = -«„ i J dE f(E - fi x „, a ) w -n a -n x ^ a (4) 

with k g = 1 (—1) and X a = S (D) for <j =^ (\). In the last passage we assumed ksT <C Eg, ^x a ,a and extended 
the integral to infinity. Note that \i is measured from the valence band edge of the TI region [Fig[2fb)]. The helical 

currents satisfy 1^ oc —f^s^ anc ^ ^[°^ ^ M^M- 

The device in Fig. [TJa) realizes the concept of a TI field effect transistor (FET), where charge transport in the ID 
helical channel can be turned on and off by changing the topological phase of a part of the system. When the device 
is ON, the source-drain current in the ID helical channel is 

AC) _ T (C) T {C) _ e 

which gives Iq ^0.1 jiA for [is = 1 meV and /i£> = 3 meV (we assume spin-balanced leads with ns,cr — l^s and 
Md,ct = Md)- This value does not significantly depend on temperature as long as /c#T < Eg, /^s, Eg — [id> 
We estimate the charge currents in the 2D bulk bands as 

J ?' ±)= ^5/* k ^dk~ lf( E -»D)-f(E-ti s )} , (6) 

where X = R or L, ± refers to conduction and valence bands, the factor 2 accounts for the spin degeneracy, Lc is 
the length of the channel and R the typical size of the contacts, a ~ * s ^ ne e ^ ec ^ 1Ye fraction of the electrons 

that has the right initial angle to propagate between the two contacts. 

We switch the device between the OFF state characterized by Vl = Vr = and the ON state with Vl = and 
Vr = 30 meV, compatible with the situation described in Fig. [2j When the FET is OFF the helical channel is absent, 
= 0, and therefore /off = ^Zx ± Iqf ^ • ^ ne ON current is instead /on = 1^ + ± Iqf ^ • ^ n Fig. [3j we 
show the current ratio /off //on as a function of temperature for different lateral size R of the contacts. It is clear 
that due to the narrow gap of the material small OFF/ON current ratios are only obtained at low temperature. We 
note however that the bulk bands responsible for the OFF current have a 2D character, while both the contacts and 
the helical channel are quasi- ID. Therefore the current ratio reduces for small R (or for long channels Lc) where 
transport in the ID edge states becomes more efficient than that in the 2D bulk bands. 

A peculiarity of the topological FET described here is that a certain amount of disorder would be favorable for 
reducing the OFF/ON current ratio. This stems from the topological protection against non-magnetic disorder of the 
helical channel, which would be little affected by the presence of inhomogeneities, while transport in the bulk bands 
(and therefore /off) would be suppressed . 
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FIG. 3: (Color online) /off/Zon current ratio as a function of temperature of the topological FET in Fig. [T] The channel 
length is L c — 1 /im, while we consider contacts of lateral size R — 0.2, 0.5, and 1 /im. 



The spin current in the helical channel of Fig. [T] is given by 

While the charge current Eq. (j5j) is evidently zero for the condition fig^ = l*>D,i, a net spin current would flow. This 
is a persistent pure equilibrium spin current, which does not lead to spin accumulation [l9[, because of the presence 
of helical edge states (even though they are gapped out) at the interface with the vacuum of the DQW. Indeed a 
conduction gap in the lateral helical edge modes is identical to a perfect backscattering process plus a spin-flip. Such 
backscattering processes lead to the following additional spin currents incoming on source and drain: 

I ( s S) = i^, 1^ = --^ D , t , (8) 

and, as expected, in equilibrium (fJLs,a = ^D,a) the total spin accumulation rate on the source (drain) due to helical 

edge modes 1^ — (Is D ^ + Is^) vanishes. 

As already noted in Ref. LL9|, |20|, spin accumulation exploiting the TI edge modes requires an applied source-drain 
bias. Let us consider the setup shown in Fig. |4j realized with an electrically tunable TI material. We assume C\ and 
C2 to be mesoscopic metallic islands (Mis), while source and drain are macroscopic reservoirs with fig^ = /is,i = l^s 
and fiD^ — {Id, i = Vd- The system is described by three characteristic times: a spin-conserving relaxation time 
te, a current injection time tj, and a spin relaxation time r s |21|. We express the particle number of the a spin 
species as ns^ = N lis, a (and similarly no, a = N /j>d,ct), with N = ^j^- the constant 2D density of states and R the 

typical dimension of the Mis. The current injection time is then defined as tj = hN = ^f-, which for rn = 0.1m e 
and R ~ 1 /im gives rj ~ 1 ns. For metallic islands of micrometer size generally te <C tj and therefore we will 
assume intra-species relaxation to be instantaneous on the time scales of 77. r s is in the range 10 -9 — 10 -6 s [2l|. If 
te <C ti <C t s the system is said to be in the slow spin relaxation limit. We set up a minimal model, which is similar 
in spirit to what was used in Refs. |2l|, [22|, where the equilibration dynamics is described by master equations, and 
neglect for simplicity charging effects on the Mis. Hence, the charge and spin dynamics in the C2 MI is governed by 
the following rate equation 

T~I Ac 2 ,a = MX CT -MC 2 ,cr [MC 2 ,o- — MC 2 ,-o-] , (9) 

while the equation for C\ is obtained by replacing /ic 2 ,a ^ ^d,-a- For spin-balanced leads with us = = ^ 

(for simplicity energy is now measured with respect to midgap) and the initially spin balanced Mis at /ic l50 -(t = 0) = 

Mc 2 ,<t(£ = 0) = /io, Eq. (0 has the solution /ic 2 ,<r(£) = ^aV + e~^" ^/io — K a Ve~~^j and /jLc 1 ,a(t) — Mc 2 ,-o-W> with 
V = Sfi 2tJ+t ' ^ n tne s ^ ow s P^ n taxation limit, /ic 2 ,o- —> k g 5fi exponentially with characteristic time tj. Opposite 
spin species accumulate on the two Mis with a spin-polarization factor P = \V/Eo\, assuming the conduction band 
edge of the Mis to be Eq < — Therefore C\ and C2 act like two opposite poles of a spin-battery. When a normal 
conducting circuit closes on such poles a pure spin current is supplied. Note however that dissipation occurs in the 
spin battery [Fig. 0] because of a finite source-drain voltage and a finite charge current. 



5 




FIG. 4: (Color online) Sketch of a spin battery realized with the exploitation of helical channels. Spin up (dashed line) and 
spin down (full line) spin filtered edge modes are shown at the boundary. 

The present configuration also realizes a set up where the dynamics of the spin densities on the MI can be studied 
by all-electrical means. In fact, by switching OFF the helical channel, once that the current is in a steady state, we 
store the spin densities on the two Mis, which then decays due to spin relaxation mechanisms. In order to measure 
the residual spin polarization fic ly a = —^c 2 ,o- = K o 5fi(t) at the delay time t, we set [is = /^d = and then turn ON 
again the helical channels. The initial value of the source-drain charge current I — — ^ 2Sja(t) allows to extract the 
residual spin polarization in the Mis. The time resolution is limited by the time it takes to switch the gate polarization 
(which does not involve any charging effects in the bulk). 

We neglected that in HgTe spin is not a properly conserved quantity but rather the total angular momentum is, 
see for example the discussion in Ref. [HI. However, it can be shown that helical edge modes are spin polarized [s[ 
(although spin is not perfectly aligned with z) and Kramers partners are characterized by opposite spin. Moreover 
the polarization axis is fixed as long as Rashba terms are negligible (23[. We also assumed that edge modes always 
ideally sink into the electric contacts, however due to a non perfect overlap of the edge mode profile with source and 
drain, they could partially avoid the contacts complicating the physical picture. For a quantitative purpose all these 
aspects should be taken into account together with the charging effects on the helical channel and, when considered, 
on the mesoscopic metallic islands. 

Many critical issues have still to be better understood before TIs can become a realistic solution for electronic 
devices. One of the limiting factors of HgTe DQWs comes from the relatively large extent (A > 100 nm) of the edge 
states, which is due to the small band gap obtainable in that system. It is therefore a crucial challenge to search 
for materials featuring larger TI gaps and well- localized edge states, like ultra-thin Bi films 0, which has been 
proposed having A < 10 nm, and devise a way to electrically tune the TI phase. Promising candidate materials include 
graphene doped with heavy adatoms [26| , silicene [27], [HI , and 2D germanium [27j • 
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